Numerical Simulation of Particle Flow in a Sand Trap 
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Sand traps are used to measure Aeolian flux. Since they modify the surrounding wind velocity 
field their gauging represents an important challenge. We use numerical simulations under the 
assumption of homogeneous turbulence based on FLUENT to systematically study the flow field 
and trapping efficiency of one of the most common devices based on a hollow cylinder with two slits. 
In particular, we investigate the dependence on the wind speed, the Stokes number, the permeability 
of the membrane on the slit and the saltation height. 

PACS numbers: 45.70.Mg,47.55.Kf,47.27.-i,83.80.Hj 



I. INTRODUCTION 



Dune motion, sand encroachment and desertification 
are based on sand transport by wind. Saltation is the ba- 
sic mechanism of Aeolian sand flux [11 |3l H] . The simplest 
and most common devices to measure it in the field are so 
called sand traps. They consist of cavities into which the 
sand bearing wind enters, drops the sand inside and then 
leaves again without the sand. In that way they accumu- 
late inside the sand that would have crossed them during 
a given time, giving a measure for the flux through their 
cross section. The difficulty to quantitative gauge them 
arises in properly estimating this cross section because 
the trap being an extended fixed object modifies consid- 
erably the wind velocity field in its surrounding. The 
issue boils down to understand how much wind actually 
enters into the cavity and therefore strongly depends on 
the geometrical shape of the trap. 

Among the many different trap designs that have been 
used in the past one of the simplest and most popular 
is a hollow cylinder with two slits, one open and the 
other covered by a membrane that is impermeable to 
the grains [1 see Fig. [l] It has been extensively imple- 
mented in field studies along the Northern coast of Brazil 
with much success. We will therefore study this particu- 
lar device in detail by numerically solving the turbulent 
flow around an object having the corresponding geome- 
try. We will calculate the trajectories of grains released 
at different positions in the area in front of the trap to 
assess if they are captured or not and do this for various 
Reynolds and Stokes numbers as well as different mem- 
brane permeabilities. This paper is organized as follows. 
In Section 2, we introduce the model and briefly summa- 
rize the algorithm used in our simulation. In Section 3 we 
present the qualitative and quantitative results. We con- 
clude in Section 4 by discussing the importance of these 
results in some practical applications. 




FIG. 1: Cylindrical sand trap in use on a field dune. 



II. MODEL FORMULATION 

As shown in Fig.|2]the sand trap that we study consists 
of a hollow cylinder of diameter D and height h closed 
on the top and the bottom, and having two vertical slits, 
one at the front and another at the back side. For sim- 
plicity we will neglect the width of the cylinder walls. 
A membrane covers the slit of the back side in order to 
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FIG. 2: The geometry of the sand trap investigated here. D 
is the diameter and h the height. The widths of the vertical 
shts at the front and back sides are Lin and Lout, respectively. 
Two possible grain trajectories are sketched. 



capture particles with diameters larger than 0.01 mm. 
The slits at the front and back sides have widths Lin and 
Lont, respectively. We chose a height h = Im, a diame- 
ter D = 0.05m and slit widths Lin = 0.01m on the front 
side and Lout = 0.02m on the back side in agreement 
with typical dimensions of real devices. The membrane 
has a thickness of Am = 0.001m. Its permeability a can 
be modified and will be one of the control parameters of 
our simulation. In order to simulate the conditions of a 
realistic wind field we consider a box of size 1.5xl.5xlm 
with moving boundaries on top, left and right, i.e. in all 
directions transverse to the wind direction except for the 
bottom. On the fixed walls, i.e. the trap surface and the 
bottom we impose non-slip boundary conditions, which 
means that the velocity is equal to zero at the interfaces 
between solid and fluid. In order to reduce the size ef- 



fects the simulation box was chosen to be one order of 
magnitude larger than diameter D of the trap. 

We assume that atmospheric air is an incompressible 
and Newtonian fluid having a viscosity of /i = 1.7895 x 
10~^kg m~^ and a density of p = 1.225 kg/m^. The 
Reynolds number is defined as Re = pVLx/fi, where V 
is the average velocity and Lx the linear size of the box 
in X direction. This gives in our case Re ~ 40000 mean- 
ing that we can assume to be in a fully developed ho- 
mogeneous turbulent state. The standard k — e model 
is therefore an adequate way to solve the correspond- 
ing Reynolds- averaged Navier-Stokes equations of mo- 
tion. Their numerical solution can be achieved by dis- 
cretizing the velocity and pressure fields and using a con- 
trol volume finite-difference technique [5 , 6 . The conver- 
gence criteria of this scheme can be defined in terms of 
residuals, which measure up to which degree the conser- 
vation laws are satisfied. In our simulation we consider 
that convergence is achieved when each of the normalized 
residuals is smaller than 10~^. 

A fully turbulent atmospheric boundary layer over a 
flat surface shows a logarithmic increase of the velocity 
v{z) with the distance z from the surface. Therefore, we 
impose a logarithmic velocity profile on the inlet, i.e. in 
the X = plane, of the form 



where = 0.36m/5 denotes the shear velocity, zq — 
1.0 X 10~^m the so called roughness length and n = 0.4 
is the von Karman constant [7 . 

We model the fabric covering the opening at the back- 
side of the trap as a special type of boundary condition 
mimicking a porous medium where the velocity /pressure 
drop characteristics are known. If Am is the thickness of 
this membrane then the pressure drop is defined accord- 
ing to Darcy's law as: 

Ap = --vAm, (2) 
a 

where a is the permeability of the medium and v is the 
velocity normal to the membrane. In order to understand 
the effect of the membrane on the wind velocity field, we 
perform simulations using different membrane permeabil- 
ities, namely a = 1.0 x 10"^, 1.0 x 10"^, 1.0 x 10"^ and 
1.0 X 10-^^ 

For simplicity we will only consider spherical particles. 
We also assume that the density of flying particles is so 
low that we can neglect collisions between them. For 
the same reason we will also neglect the momentum loss 
exerted by the particles on the fluid. Consequently the 
wind velocity and pressure fields can be calculated with- 
out knowing the particle positions and one can then ob- 
tain the trajectory of each particle by just integrating its 
equation of motion: 
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FIG. 3: The pictures show the velocity profile in and around the sand trap. The images at right, show two different velocity 
profiles in a two dimensional cut (plane) through the sand trap at height z = 0.2 and z — 0.6. The wind velocity profile input 
at the surface y — comes form Eq. ([T]). The colors indicate the magnitude of the velocity varying from blue (low value) to 
red (high value). The picture corresponds to permeabilities a = 1.0 x 10~^. 



where rup and Up 



are the mass and the velocity of a 
particle of diameter dp and ^ Fp is the total force acting 
on this particle. Let us assume that drag and gravity are 
the only relevant forces and that the particles do not 
interact with each other. Then the equation of motion 



for one particle can be rewritten as 



dwp 



{pp - P) 
Pp 



(4) 



where g is gravity acceleration and pp = 2650 kg m ^ 
is a typical value for the density of sand particles. The 
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term Fd(u — Up) in Eq. Q represents the drag force per 
unit particle-mass where 



Fd 



18/i C^Rcp 



ppdp 



24 



(5) 



Here Rcp = pdp | u — Up | //i is the particle Reynolds num- 
ber and for the drag coefficient Cd use an empirical 
expression taken from Ref. [8j. The trajectories of the 
particles are calculated by numerically integrating the 
equation of motion Eq. Q. 

We characterize the effect of inertia on the air borne 
grains in the flow field through the dimensionless Stokes 
number St which is defined by St = ppdp'^V/lSLxp. For 
St ^ 1, inertia will dominate and the particles will move 
along straight lines not following the fluid, i.e. they will 
move in the air ballistically. For St <^ 1 particles behave 
as tracers and will perfectly follow the streamlines of the 
fluid. In our simulation we will vary the Stokes number 
by changing the particle diameter keeping all the other 
parameters constant. If a particle hits the outer surface 
of the sand trap it will be reflected elastically. If it hits 
the inner part of the trap it will stop moving, i.e. is 
counted as captured. 



III. RESULTS AND DISCUSSIONS 

We ffist present a study of the influence of the per- 
meability a of the membrane on the fluid velocity field 
without the presence of particles. In order to understand 
essential features of the flow field around and inside the 
sand trap, we have performed extensive simulations by 
analyzing the flow for several values of permeability a. 
To see the flow field we focus our attention on horizontal 
cuts through the trap, at a different distances z = 0.2, 
and z = 0.6 from the ground. Fig.|3]and Fig.|4]depict the 
velocity field around the sand trap for two different val- 
ues of permeability, a = 1.0 x 10~^ and a = 1.0 x 10~^^, 
respectively. Here we focus our discussion on the ve- 
locity profile, in the horizontal cuts through the trap at 
different heights z. We will ffist discuss the case of high 
permeability. In this case we can see in Fig. |3] that the 
wind velocity presents a parabolic profile in the channel 
inside the sand trap and also the existence of stagnation 
zones beside this channel. 

For low permeability, see Fig. [4j the velocity decreases 
along the center line to form a zone of low velocity close to 
the entrance of the sand trap while it increases around 
the trap. This increase of flow is due to the obstruc- 
tion caused by the trap to the incompressible fluid. For 
larger heights the zone of low velocity reaches almost 
the entire region inside of the sand trap. With a decreas- 
ing membrane permeability also the formation of shadow 
zones can be observed behind the sand trap. 

Our central issue is the study of sand transport and 
in particular the effect of the membrane permeability on 
the sand flux. Therefore we calculate particle trajecto- 
ries by extensive simulations for different Stokes numbers 
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FIG. 4: The pictures show the velocity profile in and around 
the sand trap. The images at right, show two different velocity 
profiles in a two dimensional cut through the sand trap at 
height z = 0.6 (a) and z = 0.2 (b). The wind velocity profile 
input comes from Eq. ([T]). The colors indicate the magnitude 
of the velocity varying from blue (low value) to red (high 
value). The picture correspond to permeabilities a = 1.0 x 
10-1°. 
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FIG. 5: The 2D view of typical particles trajectories for three 
different Stokes numbers, (a) St = 0.014, {b)St = 1.43 and (c) 
St = 14.34. The membrane permeability is a = 1.0 x 10~^. 
The particles concentrate in the region of high velocities for 
all values of St. 



FIG. 6: Typical particle trajectories projected on the x — 
y) surface, for three different Stokes numbers, St = 0.014 
St = 1.43 and St = 14.34 correspond to (a), (6) and (c), 
respectively. The value of the membrane permeability is a = 
1.0 X 10"^°. 



and membrane permeabilities. In order to quantify the 
capture process the sand trap efficiency r] is defined as 



(6) 



where is the total number of particles released during 
the interval Ax and (j) is the number of particles among 
these that have been captured. We will discuss in the 
following the efficiency 77 as function of membrane per- 
meability and Stokes number. 

Initially the grains are placed at = 0, a given height z 
and randomly within the interval, —0.01m < x < 0.01m. 



For a fixed value of St and a, we release up to 1000 
particles to determine the efficiency of the sand trap. 
We show in Figs. [5] and |6] the trajectories of 7 parti- 
cles for two different permeabilities, a = 1.0 x 10~^ and 
Q/ = 1.0 X 10~^^, respectively, and different Stokes num- 
bers. For both figures the Stokes numbers are from top 
to bottom: St = 0.014, St = 1.43 and St = 14.34. In the 
case of high permeability, as shown in Fig. |5j the grains 
concentrate in the regions of large velocities. Inside of 
the trap the particle trajectories do spread apart more 
with decreasing Stokes number. This effect is more pro- 
nounced at lower permeability. For high Stokes numbers 
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FIG. 7: Semi-log plot of the efficiency 77 of the sand trap 
as function of the Stokes number. The symbols correspond 
to different values of permeabilities a, namely, (circles) a — 
1.0 X 10"^ (squares) a = 1.0 x 10"^ (stars) a = 1.0 x 10"^ 
and (triangles) a = 1.0 x 10""^°. 



some particles collide with the trap surface revers their 
direction and leave the box simulation. For intermediate 
Stokes numbers these particles also collide but are then 
deviated by the wind flow. This effect is independent of 
the permeability. 

The dependence of the sand trap efficiency r] on the 
Stokes number is shown in Fig. [7| for different values of a. 
We performed simulations for membrane permeabilities 
a = 1.0 X 10-^ 1.0 X 10-^ 1.0 X 10-^ and 1.0 x 10"^° at 
many different heights with an inlet velocity given by 
Eq.[ll The value of 77 was obtained from the average over 
different heights where at each one we released 1000 
particles. We see that r] presents two different regimes as 
function of the permeability a. Let us first discuss the 
case of low permeabilities. In the limit of small Stokes 
numbers the efficiency of the sand trap is small and re- 
mains essentially constant close to zero. Since St the 
particles can be considered as tracers that nearly exactly 
follow the streamlines of the flow, avoiding trapping be- 
cause practically no streamline crosses the interior of the 
trap as consequence of the low membrane permeability. 
Indeed we observe in these cases stagnation zones inside 
the sand trap. Above St ^ 0.05 the efficiency rj increases 
as function of the Stokes number. For high values of St 
the efficiency reaches a saturation value close to unity. In 
this limiting case the particles move ballistically towards 
the sand trap and those that have been released within 
the range —0.005m < x < 0.005m are captured. 

For high permeabilities a, the efficiency 77 presents an 
unexpected behavior in the region of low St^ as shown in 
Fig. [Tj For large Stokes numbers the efficiency 77 remains 
constant until St ~ 10.0. Below this value surprisingly 




FIG. 8: 2D view of typical particle trajectories for high per- 
meability a — 1.0 X 10~^ and low Stokes number St — 0.14 at 
height z — 0.2. The distance 5 is the maximum distance from 
the symmetry line for which the particles bend towards the 
entrance of the sand trap. The stagnation zones are shown in 
more detail in the inset revealing how the particle trajectories 
can bend towards the symmetry line of the sand trap. 



we detect a small minimum followed by an increase to a 
value that can be above one. We can explain this behav- 
ior observing the particle trajectories in the region close 
to the entrance of the sand trap as shown in Fig. |8] 

For small membrane permeability the sand trap be- 
haves like a solid cylinder. In this case a stagnation re- 
gion appears in front of the trap around the symmetry 
line. As the permeability a increases some air can cross 
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FIG. 9: Semi- log plot of the efficiency 77 as function of the 
membrane permeability a, for fixed Stokes number St = 
1.0 X 10~^. The doted line is included as guide to eye and 
corresponds to 77 = 1. 



the sand trap and two stagnation zones appear beside the 
entrance and outside of the trap. Each of them has a sep- 
aration point and all particle trajectories between these 
two separation points bend towards the center and enter 
the trap as confirmed in the simulation for low Stokes 
number as shown in Fig. [Sj Therefore the number of 
captured particles and the efficiency of the sand trap can 
increase beyond the ballistic case. 



For a fixed value of St and different permeabilities a, 
we search the maximum distance from the symmetry line 
5 at which the particles bend towards the entrance of the 
sand trap. This S can be used as a measure for the in- 
crease in the sand trap efficiency. In Fig. [9j we show in a 
semi- log plot the efficiency 77, which has been calculated 
from the parameter 6, against the membrane permeabil- 
ity a. Clearly, the curves display a strong change starting 



at the permeability a = 1.0 x 10~^. This is the smallest 
permeability which still affects the efficiency of the sand 
trap. 



IV. CONCLUSIONS 

In this paper, we have studied numerically the behav- 
ior of a sand trap frequently used to measure aeolian sand 
transport in the field. We solved the turbulent wind ve- 
locity field in the presence of the sand trap and investi- 
gated the effect of the membrane permeability. We stud- 
ied quantitatively the particle trajectories carried by the 
fluid for different membrane permeabilities. 

We have shown how the efficiency of the sand trap 
depends on the Stokes number and gave some insight 
about the effect of the membrane permeability on the 
capturing process. 

As previously observed the membrane permeability 
strongly influences the capture process at low Stokes 
numbers. The sand trap efficiency rj exhibits a surpris- 
ing increase above one, for large permeabilities and small 
Stokes numbers. Let us point out that only at higher 
St values, when the trajectories of the particles are com- 
pletely ballistic, the efficiency r] becomes independent on 
the membrane permeability. This is a useful result since 
for natural sand, for example on dunes, the grain size is 
around 200/im which corresponds to a rather high Stokes 
number. Our results therefore confirm that this type of 
sand trap is adequate to measure aeolian sand flux. The 
area of the slit at the entrance essentially equals the cross 
section over which the sand flux is measured indepen- 
dently on the details of the membrane on the back side. 
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